Perfect Actions for Scalar Theories * 

W. Bietenholz, HLRZ c/o Forschungszentrum Jiilich, 52425 Jiilich, Germany 



We construct an optimally local perfect lattice action for free scalars of arbitrary mass, and truncate its 
couplings to a unit hypercube. Spectral and thermodynamic properties of this "hypercube scalar" are drastically 
improved compared to the standard action. We also discuss new variants of perfect actions, using anisotropic 
or triangular lattices, or applying new types of RGTs. Finally we add a A(/!>^ term and address perfect lattice 
perturbation theory. We report on a lattice action for the anharmonic oscillator, which is perfect to 0(A). 



Many examples have shown that the perfect 
action program to construct improved lattice ac- 
tions works beautifully, if it can be properly im- 
plemented. However, a convincing application to 
QCD is still outstanding. Such attempts are a 
desperate struggle for locality and questions of 
parameterization and truncation are major issues. 
Hence it is motivated to study the properties of 
perfect action carefully in simple situations. 

For free and perturbatively interacting fields, 
perfect action can be constructed conveniently by 
a technique that we call "blocking from the con- 
tinuum" [0. It corresponds to a blocking factor 
n RGT in the limit n — > cx), so that the ROT 
does not need to be iterated. For a scalar we can 
relate the continuum resp. lattice field tp, (p as 

Tl{x - y)ip{y)dy , (1) 

where x e Z'^ and H(m) = 0)1=1 ©(1/2 - \u^,\). 

If we implement this relation by a Gaussian 
RGT term - with coefficient 1/a- then we obtain 
in momentum space the perfect lattice action 

Gik) = E , o-n2 , I2 + (2) 



(fc -I- 27r/)2 + m2 



where n(fc) = 11^ ^m/^a" = 2sin(fc^/2) and 
B =] — 7r,7r]''. Remarkably, the spectrum of this 



action, E'^{k) = (fc + 2TTlf 



is exactly the 



continuum spectrum, which shows that the per- 
fect action can reproduce the full Poincare sym- 
metry in observables, even though this symmetry 
is not manifest in the action. 



(3) 



If we choose the RGT parameter a = a{'m) 
(sinh(TO) — m)/ra^, then the couplings in 

are restricted to nearest neighbors in d = 1 [§, 
and they decay exponentially and very fast in d > 
1. The question arises, if this choice is really 
optimal for locality in d = 4, as we postulated 
before for fermions . We measure the decay by 
p{i, 0, 0, 0) oc exp{— c(a)«}, and Fig. |l| shows c(a) 
for various masses. Indeed, the peaks are just at 
a(m). 



a (2) 


a (4) 




a(0) 


















^--^^^ m=4 












m=2 












— m=0 















0.2 



0.4 



0.6 



0.8 1 

a 



Figure 1. The locality parameter c(a). 

As in the case of fermions, locality becomes 
even better for increasing mass, if we use Oi{m), 
see Fig. ^. 

To make this action applicable we need to trun- 
cate the couplings to a short range. We do so by 



*Talk presented at LAT97. 



^For a finite blocking factor n and m = Q^a is replaced 
by Bin = (1 — l/'T''^)/6, in agreement with B. 
^For heavy fermions, on the other hand, we noticed that 
locality can be improved slightly beyond the Id formula. 
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Figure 2. The locality measured by c(m) in 
p{i, i, i, i) oc exp(— c(m)i), for a{m) and for a = 0. 

imposing periodic boundary conditions over 3 lat- 
tice spacings, and use the resulting couplings then 
in any volume Q| . In contrast to a truncation in 
c-space, p{r) vanishes smoothly at the edges of 
a unit hypercube, the normalization is automati- 
cally correct and the "decimation in p-space" is 
handy in perfect lattice perturbation theory. 

The resulting couplings represent "smoothly 
smeared" lattice derivatives. Alternatively there 
are various ways to apply Symanzik's program to 
cancel the standard O(a^) artifacts. Symanzik 
himself suggested to use additional couplings at 
distance 2 on the axes. These couplings are com- 
pletely different from the perfect truncated ones. 
If, however, one constructs a Symanzik improved 
hypercube scalar, then the couplings look very 
similar ||^ . The same behavior has been observed 
for staggered fermions |^]. Other Symanzik im- 
proved fermions in the literature (D234, Naik) 
are constructed along the axes again. However, I 
would rather recommend to use couplings in the 
unit hypercube, which is e.g. more promising for 
the restoration of rotational invariance. 

In Fig. H we compare the dispersion relations 
of our hypercube scalar, the standard formulation 
and the Symanzik scalar on the axes. As in the 
fermionic case |^ , the latter is good at small mo- 
menta, until it is hit by a "ghost" and then the 
solutions become complex, i.e. useless. The hy- 
percube scalar, on the other hand, behaves very 
well all the way up to the edge of the Brillouin 
zone. This behavior persists at m > 0. 

The hypercube scalar has also good thermo- 
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jure 3. The spectrum along the (110) direction 
TO = 0. 



dynamic properties. Fig. Q shows the ratio 
pressure/(temperature)'', which is 7r^/90 in the 
continuum. For a small number Nt of discrete 
points in Euclidean time, this ratio is approxi- 
mated well for the hypercube scalar, whereas the 
standard action requires a large Nt to converge. 
Other quantities like the energy density look sim- 
ilar §. 
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Figure 4. The ratio p/T^ depending on Nt- 

In thermodynamics it is fashionable to use 
anisotropic lattices. There is no additional prob- 
lem to put a perfect action on an anisotropic lat- 
tice. In the most general case, lattice spacings 
(ai, . . . , flrf), wc just have to substitute in eq. (0) 
(h/ai, . . . Jd/a-d), ki_, (2/a^) sin(/c^a^/2), 
and with a —^ a/a^ the mapping on the u axis 
is ultralocal. In the typical case one has Od < 
fls = cispatiai, and One might be afraid of time- 
like "ghosts" . One can avoid them by giving up 



3 



perfection in the temporal direction, 

.^^-_^{k + 2nl/asr + kj + m^ 

We can also put the perfect action on a trian- 
gular lattice. In d = 2 we obtain a modified 11 

function, Htria{li) — 

8[fci cos ^ + cos ^ - (fci + fcz) cos 

3fcifc2(fci + /C2) 

where we refer to axes crossing under 7r/3. The 
blocking from the continuum then leads to 
3 
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Again the spectrum is perfect; we have e.g. full 
rotational invariance in the observables, although 
the lattice structure is visible in the action. 

Back to the hypercubic lattice with spacing 1: 
we can also vary the convolution function in the 
blocking from the continuum, 



Sensible candidates are the B-spline functions 
/o(s) = <5(s), /„+i(s) = r^^ljlUt)dt. The ap- 
propriate convolution function for the gauge field 
then reads fn+i{Uf,)Y{^^^, fniu^)- This set 
of functions has the correct normalization and 
a "democracy of the continuum points" , which 
all contribute with the same weight to the lat- 
tice variables. In the propagator of eq. (^, the 
power of the IT function generalizes to 2n. Deci- 
mation {n — 0) fails because the sum over I di- 
verges in d > 1, and n = 1 is what we had before. 
For higher n we can still achieve ultralocality in 
d = 1 by adding kinetic terms in a, but the lo- 
cality in d > 1 is best for n = 1. Generally, 
overlapping blocks seem to be unfavorable for lo- 
cality. We can contract /„ into [-1/2,1/2], and 
re-adjust the normalization (giving up "democ- 
racy" ) . Thus /2 turns into an "Eiffel tower func- 
tion" /2(s) = (2 - 4|s|)e(l/2 - |s|), which is in 
business in view of the induced locality H. For 



n > 2 one obtains blocking schemes which are 
difficult to relate to finite blocking factors, hence 
they are not easily compatible with the multigrid 
improvement. So I recommend the "Eiffel tower 
function" /2 as a promising alternative to /i. 

Finally we proceed to the \(f)^ theory. If we 
block perturbatively from the continuum, we en- 
counter divergent loop integrals, which can be 
regularized in the continuum. There is no di- 
vergence, however, for the anharmonic oscillator 
{d — 1). There we constructed an 0(A) perfect 
action. It involves 2- and 4-variable couplings 
oc A, which have an analytic form in momentum 
space and which we evaluated numerically in c- 
space. For a{m) they do not couple any variables 
over distance > 2 (the perfect action to 0(A") 
extends to maximal distances 2n {n > 1)). 

The performance of this action is of interest 
in view of the direct application of the perfect 
quark-gluon vertex function in QCDjT]]. We mea- 
sured the first two energy gaps |^ AEi and 
A.E2, for our action and for the standard ac- 
tion, and observed that they are clearly closer 
to the continuum results for the new action up 
to A = A/m^ ~ 0.25. This corresponds to an 
improved asymptotic scaling. However, for the 
scaling quantity AE2/AE1 the improvement is 
difficult to demonstrate, because it is restricted 
to tiny A. In both cases, the 0(A) perfect action 
is successful up to the magnitude of A, where also 
first order continuum perturbation theory - for 
the considered observable - collapses. 
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